Optimal conditions for observing Josephson oscillations in a double-well Bose-gas 
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The Josephson oscillations between condensates in a double-well trap are known theoretically to 
be strongly effected by the mean field interaction in dilute atomic gases. The most important effect is 
that the amplitude of oscillation in the relative population of the two wells is greatly suppressed due 
to the mean field interaction, which can make it difficult to observe the Josephson effect. Starting 
from the work of Raghavan, Smerzi, Fantoni, and Shenoy, we calculate the maximum amplitude 
of oscillation in the relative population as a function of various physical parameters, such as the 
trap aspect ratio, the Gaussian barrier height and width, and the total number of atoms in the 
condensate. We also compare results for 23 Na and 87 Rb. Our main new result is that the maximum 
amplitude of oscillation depends strongly on the aspect ratio of the harmonic trap and can be 
maximized in a "pancake" trap, as used in the experiment of Anderson and Kasevich. 

PACS numbers: 03.75.Fi. 05.30.Jp 32.80.Pj 



I. INTRODUCTION 

In this paper we consider a Bose-Einstein condensate 
(BEC) at zero temperature trapped in a double-well po- 
tential, which is created by superimposing a harmonic 
trap and a Gaussian barrier. With a finite barrier height 
that is greater than the condensate chemical potential, 
the Gross-Pitaevskii dynamics of the condensate can be 
well described as a superposition of right and left lo- 
calized condensates, the so-called two-mode ansatz jlj. 
The weak link allows Josephson tunneling of condensate 
atoms between wells when there is a chemical potential 
difference (or a difference in phases) between the wells. 

A double well system was implemented in the exper- 
iment described in Ref. where the Gaussian barrier 
was created by focusing a blue-detuned far-off-resonant 
light sheet into the center of a magnetic trap. The trap 
had a "cigar" geometry, that is A ^ 1, where we have 
defined the aspect ratio of the radial and axial trap fre- 
quencies A = ujp/uj z , which turns out to play an im- 
portant role in our calculations. In the experiment of 
Ref. 0]j an infinite barrier along the z-axis was used 
to separate two independent condensates of about five 
million 23 Na atoms and then removed to allow the con- 
densates to interfere. In this paper we focus on Joseph- 
son tunneling through a finite barrier, which was not 
reported in Ref. g. To date, there have been no reports 
of experimental observations of Josephson tunneling in a 
double-well trap. Josephson tunneling of a condensate in 
a one dimensional optical lattice was reported in Ref. || . 
Each well in the array held a condensate of about one 
thousand 87 Rb atoms. The chemical potential offset due 
to gravity caused the condensates to tunnel out of the 
wells. A crucial characteristic of the experiment is that 
the wells had an extremely "pancake" geometry (A -C 1) 
with A w 0.002 0,0], in contrast to the "cigar" geometry 
of the trap used in Ref. II . A related system of two hy- 
perfine states coupled by applied electromagnetic fields 
was reported in Refs. 0, Ifj, |7) for the case of strong cou- 



pling. Tunneling in a three-component condensate was 
reported in Ref. ||. 

The double-well tunneling system for a Bose conden- 
sate has been treated extensively in the recent theoretical 
literature [h], [ll], |l2[| . Nonlinear Josephson equa- 

tions describing the oscillations of the relative population 
7]{t) = (Ni(t)-N 2 {t))/N c and phase cj>(t) = <j>x (t) - ^gft) 
between the wells have been thoroughly explored |l3) , 
and this analysis will be the basis of our paper. Here, 
the indices 1 and 2 label the left and right wells and 
N c = N\ + N2 is the fixed total condensate population, 
so that — 1 < r)(t) < 1. A crucial property of this sys- 
tem is that, due to the nonlinear terms in the Josephson 
equations arising from the mean field interaction, there is 
a maximum amplitude of oscillation ?7 max < 1 less than 
unity. If the initial population difference ry(0) > ?7max, 
then very little tunneling occurs, a phenomenon referred 
to as macroscopic quantum "self-trapping" The cal- 
culations of Ref. (12) were carried out using realistic pa- 
rameters corresponding to the experiment of Ref. Q , but 
used a finite barrier height to allow tunneling to occur. 
An extremely small value of r7 max "C 0.01 was found, in- 
dicating that Josephson oscillations in that system would 
be very difficult to observe due to experimental error in 
measuring 77 (t) (which we assume to be of the order of a 
few percent). 

Because of the importance of having clear experimental 
demonstrations of quantum phase oscillations, we feel it 
is useful to explore a wider range of physical parameters 
in order to find optimal conditions for observing Joseph- 
son oscillations in a double-well. This is the main purpose 
of the present paper. The key quantities that determine 
7?max are the condensate density, the s-wave scattering 
length a, the atomic mass m, and the coupling strength 
between wells. In the absence of interactions (a — 0), the 
problem reduces to a simple linear two-level system, so 
that ?7 m ax = 1- The effect of the mean field interaction 
can be reduced by reducing the density, scattering length, 
or the mass. In this paper, we do not vary m or a, but 
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rather, we consider two types of atoms, 23 Na and 87 Rb, 
each with a fixed mass and scattering length. We calcu- 
late r? max for a wide range of trap parameters by varying 
the trap geometry determined by lo z and uj Pi the bar- 
rier geometry, and the number of condensate atoms N c . 
Wc find 77 m ax is optimized by making the trap more "pan- 
cake" shaped (A <C 1), which has the effect of maximizing 
the weak link area, while maintaining a sufficiently low 
condensate density. This result is consistent with the fact 
that Josephson tunneling was clearly observed in Ref. || 
for an array of wells possessing an extreme "pancake" 
geometry. This behavior is in contrast to that for a trap 
with a "cigar" geometry, which minimizes the weak link 
area and compresses the condensate in two dimensions 
rather than one, resulting in a greater density. 



II. DERIVATION OF MODEL 

At zero temperature, the condensate dynamics is well 
described by the Gross-Pitaevskii (GP) equation 

ih ^at^ = ^° + - 9nc(r ' ^ * (r ' tl (1) 

where n c (r,t) = |$(r,i)| 2 is the condensate density and 
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#o = -^V 2 + C/ ext (r). 
Am 



(2) 



Here, the external potential J7 xt(r) = Uh{y) + Ub(y) is 
created by superimposing a harmonic trap 

U H (p,z) = m^ 2 (AV + z 2 )/2, (3) 

and a Gaussian barrier along the z axis. 

U B (z) = U exp(~z 2 /2a 2 ). (4) 

In our calculations, we consider a range of values for the 
trap frequencies u> z and lo p , the barrier height Uq and 
width a, and the condensate population N c . 

A. Two-mode ansatz 

For a barrier height that is sufficiently larger than the 
condensate chemical potential 



(5) 



there exist two nearly degenerate stationary solutions of 
the GP equation 



[H +.gn c (r)]$ Q (r) = ^ Q $ Q (r), 



(6) 



where n c (r) = iV c |$ Q (r)| 2 . The index a = {S, A} indi- 
cates the reflection symmetry along the z-axis, symmet- 
ric or antisymmetric, respectively. These solutions sat- 
isfy the orthonormality condition J dr<f> a (r)&p(r) = 5 a p. 



Left and right localized states can be formed from these 
eigenstates $i(r) = [$s(r) + $^(r)] /\/2 and $2(1") = 
[<&s(r) — ^(r)] /\/2, respectively, which are also or- 
thonormal. 

In the two-mode approximation, the condensate dy- 
namics is described by the following ansatz solution j]J 



$(r,t) = ^i(t)$ 1 (r)+^ a (t)$ 3 (r) 



(7) 



Here, the complex coefficients ipi(t) — y/NjJt) exp[i(j>i(t)] 
are spatially uniform and contain all of the time depen- 
dence, while the two states $i(r) and "^(r) are localized 
in the left and right wells, respectively, and contain all of 
the position dependence. In a more accurate theory, one 
could account for the slow time evolution of the states 
$i(r) due to the mean field interaction since the popula- 
tion Ni(t) in each well evolves in time. If the amplitude 
of oscillation is relatively small, this can be shown to have 
a negligible effect fl pi. 



B. Gaussian variational solution 

Since we are exploring the parameter space of 
{uj z , uj p , Uo, er, N c } and must solve (||) many times, it is 
useful to construct a simplified model solution of the left 
and right states $^ (r) . Several earlier studies have shown 
that a Gaussian variational solution gives an accurate 
description of the condensate E^, [f(], 0. We therefore 
employ this procedure by taking the following ansatz for 

*«(r) 



■Mr) = C a e 



-Bl(z+z a ) 2 



± e 



-Bl(z-z a ) 2 



Here, C a is a normalization constant given by 



C a — 



7r3/2[i ±ex p(-2 J B 2 z 2 )] 



(8) 



(9) 



The three variational coefficients {A a , B a , z a } must be 
determined by minimizing the Gross-Pitaevskii energy 
functional 



dr 



^|V$| 2 + [/ cxt |<E>| 2 + f|<i>|' 



(10) 



with respect to variations in {A a , B a , z a }. We note that 
these coefficients will be different for the symmetric and 
anti-symmetric states. The Gaussian form of our varia- 
tional ansatz permits us to evaluate the integral in ( JTg| ) 
analytically. However, due to the nonlinear terms in the 
resulting expression, we carry out the minimization of 
E(Q) numerically. In Section III, we compare the Gaus- 
sian variational solution with the full numerical solution 
of ([}]) and find qualitative agreement over a wide range 
of parameters. 
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C. Equations for the relative phase and population 

We substitute the two-mode ansatz (J?]) into the GP 
equation ([!]), multiply by $i(r) and integrate over posi- 
tion to obtain equations of motion for the complex coef- 
ficients ipi(t) |J 

ihij>i(t) = [Ek + UiNi{t)]il>i(t) 

- [K + Ui,iaNi(t) + U 2 ,i2N 2 (t)]^(t), (11) 

where we have used the orthonormality property. Here, 
Ei is the bare zero-point energy in well i, 

E i= f dr$i(r)H <Z>i(r) (12) 

and UiNiit) is the mean- field interaction energy with 

Ui = g J dr$?(r). (13) 



The term K in ( |l 1| ) is the coupling energy between 
condensates due to the finite probability of tunneling 
through the barrier, 



K = 



rfr$ 1 (r)F ^2(r) 



(14) 



In (|llj), we also include the mean-field contribution to 
the coupling, given by [£/i,i2iVi(t) + U 2j x 2 N 2 ii)], where 



Ui, 12 = -g J dr$ 4 2 (r)$!(r)$ 2 (r). 



(15) 



As shown recently in Ref. |18| , these factors (|15|) give rise 
to a time-dependent coupling term, but were neglected in 
the calculation of Ref. jy. The integrals in (|l2])-(|l5|) can 
be carried out analytically using the Gaussian variational 
solution (||). 

It is straightforward to obtain the equations of motion 
for the phases and populations of each condensate by 
substituting ^(t) = yjN i {t)e i ' t ' i< ^ into (O), 



Ni(t) 



1 

h 

+ K 



Ui >12 Ni(t) + U2,i2N 2 (t) 
[NiWNjt)} 1 ' 2 sin[&(t) - <(>&)] (16) 
h&(t) = -Ek- UiNiit) + \u hX2 Nx{t) + U 2>12 N 2 {t) 

'iVj-(t)\V2 



+ K 



Using (16) and (|17|), we can derive two equations for 
the relative phase 0(f) = 4>i(t) — 4> 2 it) and the relative 
population n(t) = [Nit) - N 2 it)]/N c § 



f)it) = wjit)^l-r?it)smct>ii), 

(18) 

m = -Au J -u c rjit)-. J it)^M, (19) 



Here we have simplified the notation by defining the 
Josephson frequency 

uj(t) = [2K + N c Ui,i2(l + r?(f)) + N c U 2 ,uil - r?(f))] /H, 

(20) 

the "capacitive" frequency proportional to the mean field 
interaction 



cue = N c iUi + U 2 )/2h, 



(21) 



and the frequency due to the difference in zero-point en- 
ergies 



Alj = [£i - E 2 + iUi - U 2 )N c /2]/h. 



(22) 



For simplicity, we restrict our calculations to a sym- 
metric double- well, since the overall behavior of ?7 max 
does not depend much on a relative offset in the well 
depths. In this case, the spatially averaged quantities 
©, (If), and (§J) are equal for each well {E x = E 2 , 
Ui = U 2 , J7i,i2 = £^2,12 = U± 2 ), allowing us to make the 
simplifications Aw = and 



luj = 2{K + N c U 12 )/h, 



(23) 



which is now time-independent . Using these results, the 
equations of motion (|l8|) and (|l9| ) reduce to 



f)it) = LOj\f\^ 

0(f) = -tJcvit) 



rj 2 it) sin 0(f), 

nit) 



V^-v 2 it) 



COS (frit). 



(24) 
(25) 



In Ref. Q, Raghavan et al. give an analytic solution 
of ( p4] ) and fl25| ) in terms of Jacobian elliptic functions. 
The solution depends on the initial conditions for the 
relative number and phase, 77(0) and </>(0). From the 
solution for 77(f) we can define two maximum amplitudes 
of oscillation, one for the case 77(0) = 0, the other in the 
case of a finite 77(0) 7^ 0. For equal populations initially, 
i.e. 77(0) = 0, the maximum amplitude of oscillation 77 max 
can be obtained from equation (B5a) of Ref. Q by taking 
cfriO) = 7r/2, which gives 



1 r 

A 



2(\/A 2 + 1-1) 



1/2 



(26) 



where A = loc/uj. In the non- interacting limit (g = 0), 
loc = and hence A = 0, so that ( p6| ) gives r) max = 1, 
as it must in this limit. In the strongly interacting 
limit A ^> 1, ( p6~| ) takes the approximate form 77 max ~ 
\/2(A — 1)/A -C 1. If instead the initial population dif- 
ference is nonzero, then a slightly higher amplitude 77 c 
can be obtained compared to that given in (p6|). Taking 
0(0) = in equation (4.9) of Ref. jj]], we obtain 



77 C = 2VA- 1/A. 



(27) 



If 7/(0) > Tj c , then the solution is self-trapped. We note 
that in the strongly interacting limit A ^> 1, these two 
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quantities are approximately equal and satisfy the fol- 
lowing relation r\ c ~ -\/2f?max- We work with 77 max in our 
calculations rather than r\ c because ( 26 ) is valid over the 
whole range of values for A. 

For small-amplitude oscillations, ( pi] ) and ( pq ) can be 
linearized and combined to give 



fj(t) + u.j(loc + uj)rj{t) = 0. 



This equation tells us that the system oscillates at the 
"Josephson plasma" frequency 



yJuj(0J C + LUj). 



(29) 



In the case where the mean field interaction is large com- 
pared to the coupling energy loq/^j ^ T this reduces 
to ujjp — ^Jojjujc- Due to the mean field interaction, 
the actual frequency of oscillation between the two wells 
can be large wjp > w z , even though the coupling is weak 

LOJ < U} z . 



D. Illustration of nonlinear oscillations 

It is instructive to show the time evolution of rj{t) in 
order to visualize the behavior of the nonlinear Josephson 
oscillations and to clarify the precise meaning of ry max , ry c , 
and wjp. In Fig. 1 we graph rj(t) for three different initial 
conditions 7/(0) and <M0) for the solution of the coupled 
nonlinear equations (|2J) and (|25|). We first consider the 
case of equal populations rj(0) = and an initial rela- 
tive phase of cf>(0) — 7r/2, shown by the solid line. The 
population difference r](t) oscillates with an amplitude 
of r) max . In the non-interacting limit (eg. luq — * 0), the 
problem reduces to a linear coupled two-level system and 
a complete population transfer between wells would oc- 
cur in this case, rj max — 1, at a frequency of u>j. The 
amplitude suppression introduced by the mean field in- 
teraction arises from the term —LOcTi(t) in (25), which 
acts as a time-dependent offset in the zero-point energies 
of the two wells. We also indicate the period of small- 
amplitude oscillation Tjp = 27r/wjp on the graph, which 
is valid when the amplitude is much less than 77 max (we 
have verified this in calculations not shown here). Tjp 
gives a useful estimate of the period, which, in general, 
depends on the initial conditions. 

We also consider the case of unequal initial populations 
7y(0) ^ 0. For the dashed line in Fig. 1, we have taken 
77(0) = 0.9r] c and cj)(0) = 0, where the critical amplitude 
is rj c = 0.11 for this case. If the initial population dif- 
ference is larger than rj c , then the population becomes 
"trapped", as shown for the case of 7/(0) = 0.15 by the 
dot-dashed line. It is important to realize that in the 
non-interacting limit (toe — ► 0), this self-trapping does 
not occur. In the absence of interactions, a full oscillation 
of the dot-dashed line would occur between 77 = ±0.15, 
at a frequency of ujj. We remark that for 77(0) ~ rj c , 
higher harmonics appear in the solution due to the non- 
linearity [0. We refer the reader to Ref. for a more 



0.15 



(28) S 



i \ 1 > 

i \ i \ 




100 \ 200 

T jp t (ms) 



300 



400 



FIG. 1: Amplitude-suppressed Josephson oscillations of r/(t). 
The nonlinear equations (^) and ( |25| ) for rj(t) and <j>(t) are 
solved for different initial conditions: 77(0) = 0, 4>(0) = n/2 
(solid), 7/(0) = 0.1, 0(0) = (dashed), r/(0) = 0.15, 0(0) = 
(dot-dashed). We indicate along the vertical axis the critical 
value 77 c = 0.11 and the maximum amplitude 77 max = 0.075 
obtained when 77(0) = 0. Along the horizontal axis we also 
indicate the estimated period of the Josesphson oscillations 
Tjp = 2n/ujjp. The physical parameters used correspond to 
the trap of Ref. § in Table 1, but for iV c = 1000 atoms. 



detailed discussion of this system, where a full range of 
solutions of (^) and ( p5| ) is explored. We also note that 
this same behavior is also found in Ref. [[[J] for a conden- 
sate with two internal states coupled by a weak external 
electromagnetic field. 



III. RESULTS 

In this section we calculate 77 max as given by (|2^) as 
a function of the trap frequencies lo z /2tt and w p /27r, the 
barrier height Uq and width er, and the condensate pop- 
ulation N c . To do this, we solve for the left and right 
localized states §i(r) and then calculate the quantities 
(|i^)-(|l5|), so that the frequencies ojj and uic can then be 
computed. We compare two different solutions of "^(r), a 
direct numerical solution |l9) of the stationary GP equa- 
tion (J6j) and the Gaussian variational ansatz (||). The 
variational coefficients {A ai B ai z a } are computed nu- 
merically by minimizing the energy functional E($>) (|Io| ) 
using the Nelder-Mead simplex search algorithm in the 
software package MATLAB. 

While varying the physical parameters lu z , lu p , ct, and 
N c in our calculations, we keep Uq at a fixed height rela- 
tive to the chemical potential us, which also varies. The 
bottom of the double well J7 m i n is obtained by minimizing 
U cx t along the z-axis to give 



Umin = muj z a 2 [l — ln(mwf a 2 /Uq)] . 



(30) 
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TABLE I: Table of trap parameters. In the calculation of 



Vn 



we take N c 



10 4 



1.5, and a = 1/im. For the 



scattering lengths of 23 Na and 87 Rb, we take dNa = 3 nm and 
«R.b = 5.7 nm. 



trap 


lu z /2tv (Hz) 


ujp/2-n (Hz) 


A 


atom 




Ref. J2 




19 


250 


13 


23 Na 


0.036 


Ref. § 




65 


24 


0.37 


87 Rb 


0.14 


Ref. |p 




8 


8 


1 


87 Rb 


0.26 


Ref. J3 




7200 


17 


0.002 


87 Rb 


0.85 



In our calculations, we define the ratio \ of the barrier 
height Uq to the chemical potential [is of the symmetric 
solution, each measured relative to the bottom of the well 



X 



(Up - t/min) 
(MS - E/min) ' 



(31) 



We emphasize that in all of our calculations, the coupling 
due to tunneling is weak, with 0.05 < u}j P /lu z < 0.5. 

We first calculate r7 m ax for a few typical magnetic traps 
shown in Table 1. For the "cigar" trap of Ref. j^j, tak- 
ing N c = 5 x 10 6 and a barrier width of a = 6/im, the 
maximum amplitude of oscillation is miniscule r/max = 
5.2 x 10~ 6 . Even if we greatly reduce the condensate den- 
sity by lowering the population iV c = 10 4 and increase 
the coupling by reducing the barrier width a = 1/im, 
the maximum amplitude is still very small, with r/ ma x = 
0.036. The Josephson oscillations would not be visible 
in this case if the experimental error in measuring r/ ma ,x 
were a few percent. For the "pancake" trap of Ref. ||, 
taking a modest condensate population N c = 10 4 and a 
narrow barrier a = 2/xm, we find r/ m ax = 0.054. A better 
result can be obtained by decreasing the barrier width 
cr = 1/im, which gives r/max = 0.14. Finally, we consider 
the spherical trap of Ref. , which is a much shallower 
trap than the other two, so that the density of the con- 
densate is much less. Taking N c = 10 4 and a = 2/xm, 
we find r/ ma x = 0.12, and for a = 1/im, 7/ ma x = 0.26. For 
the optimistic parameters of this shallow trap, Josephson 
oscillations should be observable. 

For comparison we also consider the experimental 
setup of Ref. ||, where a condensate was loaded into 
a one dimensional optical lattice. Each well in the ar- 
ray can be modeled as being harmonic as a first ap- 
proximation, for which one obtains lu z /2tt ~ 7200 Hz 



and u> p /2ir ~ 17 Hz 



For a rough comparison, 



we consider a double well with this geometry. Taking 
X = 1-5, a = 1/im, and N c = 10 4 87 Rb atoms, we find 
?7max = 0.85. This value increases to r/ max = 0.99 if we 
reduce the population to N c = 10 3 . 

Motivated by these results, we obtain a very crude es- 
timate of the scaling of r7 m ax with the parameters. In the 
strongly interacting limit A ^> 1, (P6|) reduces to 



Vn 




(32) 



1 0.6 




FIG. function of lu z /2tt for three values of uj p /2n. 

The three lines are obtained using the Gaussian variational 
solution (|^) and the filled circles are obtained from a full nu- 
merical solution of the GP equation (^|). In this calculation, 
the condensate population is held fixed at N c = 10 3 atoms 
of 23 Na. The relative barrier height ( fjl] ) and width as given 
in (0) are fixed at x = 1-5 and a = 1/im, respectively. 



Taking a simple WKB approach to calculate the tun- 
neling probability p] , pp[ , an approximate scaling of 
the Josephson frequency can be obtained luj/lu z ~ 

cxp[— y ma 2 {Uo — /i) /h 2 } . One can also get the approx- 
imate scaling of u>c on m and a by taking u>c ~ A*tf, 
where /xtf/^Wz = 0.5[15A 2 A f c a/zh o ] 2 ^ 5 is the chemical 
potential of a condensate in a harmonic trap in the 
Thomas- Fermi limit pl| . Inserting the definition of the 
oscillator length, Zho = yS/mWj, the expression can 
be simplified to luc/^z ~ (ma 2 N 2 \ 4 uj z /h) 1 ^ 5 . Insert- 
ing these crude estimates of loj and loc into ( |32"| ) , we find 
that r/ ma x scales like 



cxp [- \Jma 2 {U - [i)/4h 2 
(ma 2 X i N 2 cj z /h) 1 / 10 



(33) 



This indicates that the amplitude of oscillation can be 
increased by decreasing m, a, N c , a, or Uq. For a fixed 
Wjt, r/max increases as A is decreased, while for a fixed as- 
pect ratio A, r/max decreases when the trap is tightened 
by increasing lo z . Expression ([33]) only gives a very crude 
estimate of the scaling behavior. In the following discus- 
sion, we explore the dependence of r/ max on the physical 
parameters more quantitatively by calculating ljj and u>c 
numerically. 

In Figs. 2 - 5 we plot r/max (given by (p6|)) as a func- 
tion of the trap geometry (by varying lu z and w p ), N c , 
and the barrier height x an d width a. In each figure, we 
show results for the Gaussian variational solution (lines) 
and the direct numerical solution of the stationary GP 
equation (0) (filled points). The Gaussian variational so- 
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'°9io N c c (urn) 



FIG. 3: ?7max as a function of the condensate population N c 
for 23 Na and 87 Rb. The lines are obtained using the Gaus- 
sian variational solution and the filled circles ( 23 Na) and 
squares ( 87 Rb) are obtained from a full numerical solution of 
the GP equation (n) . In the inset we show how the frequency 
of oscillation lojp/lo z varies with N c . In this calculation, the 
harmonic trap frequencies are held fixed at v z — 100 Hz and 
v p = 10 Hz. The relative barrier height and width are fixed 
at x — 1-5 an d a — lM m > respectively. 



lution agrees with the full solution within a few percent 
over a broad range of parameters. The greatest discrep- 
ancy occurs in Fig. 2 as ui z decreases at fixed oj p and 
in Fig. 4 as a increases, with the two solutions differing 
by a factor of 2 at a = 4/xm for 87 Rb. In all cases, the 
Gaussian solution slightly underestimates r) max . 

Our most important result is shown in Fig. 2, which 
shows the strong dependence of ?7 max on the trap geom- 
etry. For a fixed value of uj p , ??max can be optimized by 
increasing the axial trap frequency u> z , that is, by mak- 
ing the trap more "pancake" like. In contrast, for a fixed 
value of uj z , r] ma , x decreases as uj p is increased. It is im- 
portant to realize that both quantities u>c and wj are 
effected by changing the aspect ratio of the trap. By 
making the trap more "pancake" like, one is increasing 
the area of contact between the two wells, which increases 
ljj. On the other hand, compressing the condensate in 
the z dimension also leads to a higher density, causing loc 
to increase as well. However, the ratio of these two quan- 
tities in (^2|) increases for increasing u> z . In contrast, by 
making the trap more "cigar" like, one is decreasing the 
area of contact between wells, resulting in a decrease in 
ujj as LOp increases. Furthermore, in this case one is com- 
pressing the condensate in two dimensions rather than 
one, which results in a more dramatic increase in the 
density and a corresponding increase in wc, so that the 
ratio ll>j/u>c decreases with increasing u> p . 

In Fig. 3 we vary the condensate population N c from 
10 2 to 10 5 atoms, and show r? roa x f° r both 23 Na and 87 Rb. 
We see in Figs. 2-4 that a much larger amplitude of 



FIG. 4: 7? ma x as a function of the barrier width a for Na and 
87 Rb. In this calculation, the harmonic trap frequencies are 
held fixed at v z = 100 Hz and v p — 10 Hz. The relative barrier 
height is fixed at x = 1-5 an d the population is N c = 10 4 
atoms. 



oscillation 77 max can be achieved using 23 Na compared 
to 87 Rb. This behavior is expected based on the rough 
scaling of ?7 max given in (33). We can understand why 
23 Na is favorable over 87 Rb, since the mass m of 23 Na 
is a about a factor of four smaller than that of 87 Rb, 
and the scattering length a of 23 Na is almost half that of 
87 Rb. As N c increases in Fig. 3, the effect of the mean 
field interaction becomes more pronounced causing loq to 
increase, while this has only a small effect on uij. We also 
show the dependence of wjp on the condensate number 
N c in the inset. 

We also show the dependence of ?7 max on the Gaussian 
barrier parameters by varying the relative height x ( |3l| ) 
in Fig. 4 and the width a in Fig. 5. We note that varying 
these two parameters has very little effect on ljc, which 
remains roughly constant in these plots. The behavior of 
f?max is not surprising: as the barrier height and width are 
each increased, the coupling between wells is diminished 
so that ujj decreases, resulting in a decrease in 77 max as 
indicated by @. 



IV. CONCLUSIONS 

In this paper we have focussed on an important prop- 
erty of the double-well Josephson system for a dilute 
atomic Bose-Einstein condensate: the maximum ampli- 
tude of oscillation 77 max of the population between wells 
is suppressed due to the mean field interaction [jl], 
O For typical magnetic trap parameters, we have 
shown that 7? max is very small, indicating that Josephson 
oscillations may be visible only for the very shallow trap 
used in the experiments described in . We have shown 
that much larger values of 77 max can be achieved by de- 
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FIG. 5: 77mas-.as a function of the relative barrier height \ 
defined in ( pl| ) for 23 Na and 87 Rb. In this calculation, the 
harmonic trap frequencies are held fixed at v z = 100 Hz and 
v p — 10 Hz. The barrier width is fixed at a — 1/im and the 
population is N c = 10 4 atoms. 

creasing the aspect ratio A = o; p /w z , that is, by making 
the trap more "pancake" like, in particular, for small val- 



ues of the radial frequency uj p . Our study was motivated 
by the fact that Josephson oscillations were clearly ob- 
served in the experiment reported in Ref. || , where a rel- 
atively small condensate was loaded into an array of wells 
with an extremely small aspect ratio A 0.002 Q. 
We have also shown that larger values of ry m ax can be at- 
tained by decreasing the condensate population N c , and 
decreasing the barrier height x an d width a. Further- 
more, conditions are much more favorable for 23 Na com- 
pared to 87 Rb. It should be noted, however, that present 
traps for 23 Na atoms possess the less optimal "cigar" ge- 
ometry. 

We have restricted our calculations to zero tempera- 
ture. Our study of the optimal conditions for observ- 
ing Josephson oscillations should be extended to finite 
temperatures, where dissipative effects associated with 
the thermal cloud of non-condensate atoms must be in- 
cluded jyj |22| . Very few explicit calculations of damping 
rates of Josephson oscillations in a dilute Bose gas have 
been carried out, and further study of damping in this 
system is needed. 
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